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Abstract. We give a quantum gate construction - composed entirely from 
incidents of the CNOT gate - that generalises the qubit SWAP gate to higher 
dimensions. This new construction is more regular than and is an improvement 
on the WilNOT quantum gate construction. 



1. Introduction 

It is becoming increasingly evident that much effort continues to be made into 
finding optimal quantum circuits in the sense that for the given gate library there 
is no smaller circuit that achieves the same task. A reason for this concerted effort 
is primarily due to quantum error which is called decoherence. While a quantum 
computer is predicated on the undisturbed evolution of quantum coherences, deco- 
herence represents a major and unavoidable problem for the practical realisation of 
quantum computers as it describes the error state of a quantum computer which is 
introduced through the interaction of system and environment. As with classical 
computations, the properties of experimentally realisable quantum gates influence 
the execution time of quantum algorithms. It is therefore advantageous to ensure 
a minimal use of computational resources to limit execution time and exposure to 
the quantum environment. Thus, a minimal use of computational resources serves 
to cap the total decohering time delivered by the execution of quantum gates. 

While universality of quantum gate constructions is a key concern, researchers 
now focus on searching for universal n-qubit gates that contain fewest uses of CNOT 
gates (Bergholm et al. (2004, 2005), Sedlak M and Plesch (2004), Shende et al. 
(2006)). This is because the cost of experimentally realising a CNOT gate exceeds 
the cost associated to single qubit gates (Bergholm et al. (2004)). Consequently, 
construction of quantum circuits that minimise the use of CNOT gates are impor- 
tant from the point of view of execution time of the corresponding quantum algo- 
rithm which is positively correlated with the decohering time of the circuit. Among 
universal designs for two-qubit quantum gates there are those containing six, four 
and three CNOT gates (Bergholm et al. (2004), Vatan and WiUiams (2004), Vidal 
and Dawson (2004). However, characterising the exact CNOT complexity of an 
arbitrary n-qubit operation and constructing the corresponding efficient quantum 
circuit is seen as an ambitious task, even by numerical analysis (Vidal and Dawson 
(2004)), but researchers have done considerable work optimising their constructions 
(Nielsen (2005)). In particular, Vatan and WiUiams (2004) give a quantum circuit 
construction for a general two-qubit operation that requires at most three CNOT 
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gates and fifteen one-qubit gates and show that their construction is optimal. Cru- 
cial to in securing this result is the demand that quantum circuits for the two-qubit 
SWAP gate composed using CNOT gates requires at least three such gates - to show 
the correctness of this claim, Vatan and Williams utilize the notion of entangling 
power introduced by Zanardi et al. (2000). Importantly, a scheme to realise the 
quantum SWAP gate qubits was presented by Liang and Li where it is maintained 
that experimentally realising the quantum SWAP gate is a necessary condition for 
the networkability of quantum computation (Liang and Li (2005)). 

Although universal two-qubit circuits with fewest uses of CNOT gates are known, 
the threc-qubit TofFoli gate is known to require at most six CNOT gates (Nielsen 
and Chuang (2000)) but only five CNOT gates have been shown to be necessary 
(Vidal and Dawson (2004)); see also references therein. Constructing a quantum 
circuit for a given n-qubit operation is seen as an important task for the complete 
networkability of quantum computation. Unfortunately, the efficiency of universal 
designs that accomplish this task is often known for the worst case of n-qubit 
operators. Although, the asymptotic number of CNOT gates used by Barenco's 
decomposition to implement any n-qubit operation is 0(n'^4") (Sedlak and Plesch 

(2004) ); it is believed that there are interesting operations which might require a 
polynomial number of CNOT gates. 

Devising quantum gate constructions that extend to n-qubits remains an open 
problem. Furthermore, it is believed that the study of quantum circuit minimisa- 
tion and generalisation of qubit circuit architectures to qudit circuit architectures 
will require different construction techniques than those presently known. There is 
another direction in which this problem may be taken. This is the generalisation 
from using qubit gates to using qudit gates for circuit design. To this end, a geo- 
metric approach to quantum circuit minimisation has been put forward by Nielsen 

(2005) that seeks to find the length of a minimal geodesic with respect to a suitable 
Finslcr metric. A striking feature of this approach is that once an initial position 
and velocity of the geodesic are determined the remainder of the geodesic can be 
completely evaluated by a second order differential equation. This is in contrast 
with the usual case of classical and quantum circuitry design where part of a circuit 
does not aid the complete design. 

Considering quantum network designs within the qudit setting taken with new 
approaches to circuitry design ought to be merit in itself. Therefore, a question 
that remains is whether a generalised SWAP gate for higher dimensional quantum 
systems can be completely and optimally constructed from uses of a generalised 
CNOT gate. Should this be the case then such a construction may be a cornerstone 
in the networkability of quantum computers based on qudits. 

2. Preliminaries 

Given an arbitrary finite alphabet E of cardinality d, we process quantum in- 
formation by specifying a state description of a finite dimension quantum space, 
in particular, the state description of the Hilbert space C*. While the state of an 
d-dimensional Hilbert space can be more generally expressed as a linear combina- 
tion of basis states \'4'i): we write each orthonormal basis state of the d-dimensional 
Hilbert space to correspond with an element of Z^. In this context the basis 
{|0) , |1) , . . . , Id — 1)} is referred to as the computational basis. Therefore, a state 
IV") of is given by \il>) — YltZo oii N)j where G C and Xlf=o^ Wif' = 1- A qudit 
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describes a state in the Hilbert space C , and the state space of an rt-qudit state is 
the tensor product of the basis states of the single system C^, written Ti, = (C'')'*", 
with corresponding orthonormal basis states given by \ii) |?2) (8) ■ ■ ■ (81 \in) — 
\i1i2 ■ ■ - in), where ij e Z^. The general state of a qudit in the Hilbert space Ti, is 
then written 

(2-1) 1^) = 'Y^ a;(,j,2...i„) Ki*2 • ■ • in) , 

where £ C and ^ 10^(1112. ■in) P ~ 1- The qudit representation of a quan- 

tum state provides a natural mechanism by which quantum computations can be 
implemented. That such a computation is made possible initially lies with the no- 
tion of state signature. In particular, the correspondence of quantum information 
ttfe with a computational qudit basis clement |fc) and the subsequent genesis of the 

quantum state X]fe=o '^'^ 1^) ™ ^^'^ Hilbert space C*. Such a correspondence between 
information and a Hilbert space representation is prerequisite to quantum compu- 
tation since the successful transmission of any information state is predicated on 
encoding the basis states associated with the quantum information elements rather 
than the information itself. 

Let 7i_4 and TLb be two d-dimensional Hilbert spaces with bases and , i S 
Z(j respectively. Let denote a pure state of the quantum system T-Ia- Similarly, 
let |(/))g denote a pure state of the quantum system Tig and consider an arbitrary 
unitary transformation U G U(d^) acting on Ha ® Hb- Let J/cnot denote a 
controUcd-NOT (CNOT) gate that has qudit \tp)j^ as the control qudit and as 
the target qudit; then 

(2.2) UcNOT \m) ^ (g) \n) = \m) ^ (g) \n ® m) , m,neZd 

where i G) j denotes modulo d addition. 

3. A Generalised Quantum SWAP Gate Construction 

Let us consider a set of d qudit quantum systems, the first system ^0 prepared 
in the state |eo)g, the second system Ai prepared in the state \ei)^ and so forth, 
with the final system Ad-i prepared in the state \ed-i)d-i- We ask whether or not 
it is possible to construct a network to implement a generalised SWAP gate so that 
in the output of the network the system Ao is in the state lei)^, the system Ai is 
in the state {€2)1 and so forth, until the system Ad-i is in the state \eo)d-i where 
Co, ... , Cd-i S {0, . . . ,d — 1}. To this end, we give the following. 

3.1. The Construction. Let k and I be positive integers. For non-negative inte- 
gers j, let us consider the function f{j) = (^) mod I and note that f{j) is a cyclic 
function of j. We use the periodic property of this function to study the com- 
putational problem associated with the construction of a generalised SWAP gate 
for quantum systems over dimensions d. In particular, we will relate the binomial 
summation aj = X]i=o ("'"'' mod d to a quantum network for a quantum gate, 
where such a network is constructed from the periodic application of CNOT gates 
connecting adjacent systems of the network. 

Fig. [1] illustrates such a quantum network for 4-dimensional quantum states. 
The diagram shows the sequence of CNOT gates of the network. The columns of 
the array on integers (modulo 4) describe the states of the target systems after 
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aj (mod 4) 



ei, 62,63, eo, 60 © 61, ... 

000111123013223202133121013001 
100011112301322320213312101300 
010001111230132232021331210130 
001000111123013223202133121013 
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Figure 1 . Binomial summation network construction over dimen- 
sion 4. 



the corresponding CNOT gates have been apphed, so the arrow indicates the cor- 
respondence between the system Ai being in the state |eo -I- ei)-^, where eo -I- ei is 
calculated modulo 4, and the appHcation of the first CNOT gate. The sum eo -I- ei 
may be represented as the dot product of the row vector 60, 61, 62, 63 and the corre- 
sponding column (1, 1, 0, 0)"^. Thus the array of integers modulo 4 has rows indexed 
by {0, 1, 2, 3} and columns indexed by time t = —3, —2, —1, 0, 1,2,3,... such that 
column t = As + j with j £ {0, 1, 2, 3} corresponds to system Aj. We denote the 
entries in the array by bit, where i G {0, 1, 2, 3}, and t = —3, —2, —1, 0, 1, 2, 3, ... . 
Putting at = 6i&it, t = —3, —2, —1, 0,1,2,3,..., then, for t = As + j with 

j e {0,1,2,3}, the state of the system Aj is \a4s+j)^ - Aj has been the target of 
s + {1 — Soj) CNOT gates. In general, for a network of d systems of dimension d, 
the periodic arrangement of CNOT gates with control system Aj and target sys- 
tem Aj+i for j = 0, . . . , c? — 1, where j + 1 = for j = d — 1, means that the 
sequence (at), where \ads+j)j is the state of system Aj at time t ^ ds + j with 
j G {0, . . . , (i — 1} and Aj has been the target of s -|- (1 — Soj) CNOT gates, satisfies 
the recurrence at+d = at+d-i + at for all t > —d + 1, since the CNOT gate replaces 
the state \ads+j)j oi Aj with \ad{s+i)+j)j = \ads+j + ad(s+i)+j-i) y With initial 
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States |a_d+i)i = , \a-i)d-i = |ed-i)d_i, |ao)o = |eo)o the terms of the 

sequence (at) may be written as at = '^'iZo ^i^u for sequences (bn) which satisfy 
the recurrences ^^(t+f;) = &i(t+(j-i) + ^iti « = 0, . . . , d — 1. Indeed 

d-l d-l 

at+d-i + at — ^ ei6i(t_|_(i_i) + ^ e^bit 

i=0 i=0 
d-l 

(3.1) = ^ej(6j(t+d-i) +feit) 
and 

d-l 

(3.2) ot+d = eibi{t+d)- 

4=0 

We also note that, for i = 0, . . . , d— 2, the sequence (6^4) is a translate of the sequence 
by 1 place. Indeed these sequences are all translates of one another. We 
begin by considering the solution of the recurrence relation a[t+d) ~ ot+d-i + a*- 
This turns out to be at = Y!J^o 

Lemma 3.1. (Rosen (2000)) (^) + = (^+^). 

Lemma 3.2. Let d he a positive integer. The sequence (a„) defined by a„ = 
X]r=o satisfies the recurrence relation an+d = On+d-i + On with initial 

conditions ao = ■ • • = fld-i = 1- 

Proof: Clearly the sequence (a„) as defined satisfies oq = • • • = a^-i = 1. Let I 
be a non-negative integer and let m G {0, . . . , d — 1}. Then 



(3.3) flid+jri — ^ 



i 



and 



a/d+rn+d-l = a(i+l)d+m- 1 = 

4=0 

( + 1 

i=i 

(3.4) = 1 + ^ 



''(Z + l-i)(i + m + i-l 



+ 1 - z)(i + TO + i - 1 



i=0 



if — i)d + m + i 
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Hence, by lemma ISTTI we have 



aid+m + dld+m+d-l 



l + l 

1+E 



' ^{l-i)d + m + i + l 
i + l 

{l + l- i)d + m + i 



i=l 

l+l 



E 



i=0 



{l + l-i)d + m + i 
i 



(3.5) = a, 



{l + l)d+m 

as required. 

The solution sequence (oj ) is periodic since the recurrence relation is reversible and 
the sequence must repeat as soon as d consecutive terms, of which there are only 
finitely many possibilities, are repeated. To determine the period of this sequence 
we need some combinatorial results. 

Lemma 3.3. (Rosen (2000)) ^Lo {'V) = C^fe^^) V ^ 0- 

Lemma 3.4. Let p be prime and let j > —1. Then for j ~ p — 1 {mod p), we have 
(pil^i) ^ 1 {mod p) and for j ~ 0, . . . ,p — 2 {mod p), we have (p^{) = {mod p). 

Proof: Let us write (p^^) as a quotient of factorials, and consequently, we have 
it that = . Since there is a multiple of p in the numerator 

not cancelled by a factor in the denominator except when j = p — 1 (mod p) , the 
result follows. 

Theorem 3.5. (Lu and Tsai (2000)) Let p be a prime number and let a, and k be 
any positive integers. The integer function {-Q modulo p"" for j > k has the cycle 
length p"-^'^ where e = {log^kl . 

We now determine the period of the solution sequence {oj) when d is prime. 

Theorem 3.6. Let d be a prime and consider the sequence {oj) defined by Oj = 
Ei=o C^^I^^^O "^0*^ ^- ^^en (flj) has the cycle length - i. 

Proof: The sequence (oj) satisfies the recurrence Oj+d = aj+d-i + for all j. 
Since oq = ai = • • • — a^-i = 1, it is sufficient to show that a(j2_j = l,a^2_2 = 
0,...,a^2_j^ ~ 0. This implies that ajj^_j2_i = 1 for j ~ 0,...,d—l, and so 
flj+d^-i = «j for all j > 0. Now, for j = 0, . . . , d — 2, 

-12 



V- fj + d^~d-{d~ l)i\ 
aj+d^-d = 2^ I • ) 

d-l 

(3.6) = E 



i=0 



j + i 
i 
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By lemma [331 



r mod d for < j < - 1 
[ 1 mod d for j — d — 1. 



Thus the period divides — 1. Next, we show that Pa = — 1 is the smallest cycle 
length. For i = 0, (^"'''~^'') 1 for aU j > 0. Let i be such that 1 < / < 
Then the sequence (cj) with Cj = ("' '^t^^'*') satisfies Cj = for j = 0, . . . ,di — 1 
and Cj = 1 for j = di and is periodic with period d for j > di by Theorem 13.51 
Thus, aid = / + 1 for ? = 0, . . . , d— 1. Now, any d consecutive terms of oq, . . . , aci2_i 
includes a term aid for some I and thus there cannot be d — 1 consecutive and one 
1 until the terms from d^ — d to d-^ — 1. This establishes the result. 

Since the cycle length d^ — 1 is coprime to d, note that at the completion of a cycle the 
states of the system will be cycled round. In fact, asd^ — l = d— 1 = — 1 (mod d) 
they are shifted by one position and the network implements the SWAP gate - 
system will be in the state |ei)Q, the initial state of Ai, Ai will be in the state 
62)]^, . . . , Ad-i will be in the state \eo)d-i- This can be seen by the following 
argument. Since the sequences (bu), i = 0, ... ,d — 1, have cycle length d^ — 1, 
we have (&ot, • ■ ■ , ^(d-i)t) equal to: (0, 1, 0, . . . , 0) for t = d^ — d corresponding to 
system (0, 0, 1, 0, . . . , 0) for t = d^ — d + 1 corresponding to system Ai; . . . ; 
(0, . . . , 0, 1) for t = d^ — 2 corresponding to system Ad-2', (1,0,..., 0) for t = d^ ~1 
corresponding to system Ad-i- 

Now let us consider the prime power dimension d = p™ with p prime and let us 
further consider the family of binomial summations X]i=o {^~^^ i ~^^'^')- We conjec- 
ture that the cycle length of the integer sequence (oj) for aj = X]i=o 
mod is p™-i(p2m _ see Tabled] Note that if this conjecture is true then 
since gcd(p™,p™^^(p^'" — 1)) = p™^^ ^ 1, the network does not produce a SWAP 
gate when m > 1. Although at the end of a cycle the systems are shifted round, 
they are shifted p™~^ places and the systems return to their original states after p 
applications of the network. Thus, the network provides a cyclic swap on each of 
p-m-i gj-Q^pg Qf p systems, the systems of a group have indices congruent modulo 
pin-i^ What distinguishes this network from copies of the SWAP network for 

d = p is that in this larger network the systems that are swapped do not directly 
interact through a CNOT gate. 

We turn to the problem of determining the cycle length of the sequence (aj) for 
d = p™, (to > 1). We have it that over such dimensions the sequence (aj) satisfies 
the recurrence aj_|-pm = aj_|_pm_i + aj for all j. Since ag = ai = • • ■ = apm_i = 1 
by definition, it is sufficient to show that apm-i(^p2m_i^ = 1 and apm-i(p2m_i)_j 
vanishes for j = 1, . . . ,p™ — 1. It then follows that apm-i(p2m_i)_|_j = 1 for j = 
0, . . . ,p™ — 1. Considering this problem, we seek a closed form expression for the 
integer sequence aj = J2i=o ("'~^^ , ~^^*) . We outline the approaches we have taken. 



d-l 

E 

i=0 



j + i 



(3.7) 
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Firstly, following the approach for d = we find that 



ap™-i(p2™_i)_j 



E 
E 

i=0 



mod 



mod 



(3.8) 



E 

im-l _ 

e" 

i=0 



d71 
P 



mod p™. 



We ask if there exists a combinatorial approach which illustrates that the sum 
^^^Q ~^ (^(i-p )'-P -i) mod equals 1 for j = and vanishes for j = 
1, . . . — 1. On the other hand, many combinatorial problems have been resolved 
through the theory of hypergeometric series (Petkovsek et ai). A hypergeometric 
series X]fe>o is one in which Iq = 1 and the ratio of two consecutive terms is 
a rational function of the summation index. The process of identifying a given 
hypergeometric series by writing the series in the standard pFq [■] may may aid a 
simple closed form for the series by comparing it with the library of closed forms 
associated with hypergeometric series listed in Petkovsek et ai. Interestingly, the 
sequence (aj) can be identified as a hypergeometric series of the form 



(3.9) 



-j+i 

pin : 



-3+P -1 
p'" 



-3+p'"-2 



(p" 



Unfortunately, the library of closed forms listed in Petkovsek et ai does not pro- 



Instead we consider the fol- 



vide a closed form for the sum 
lowing approach for evaluating the closed form of the binomial summation aj = 

Ei^^'C-^^T-'^Vodp™. 

Let A{z) be the power series J2j>o '^i^'' denote by [z^]^(z) the coefficient of 
in A{z)] thus [z-']^(z) = aj. Now we determine the generating function for A{z). 
A generating function is a clothesline on which we hang up a sequence of numbers 
for display (Wilf). In particular, the jth term of the sequence (oj) is the coefficient 
of z^ in the expansion of its generating function as a power series. To find the 
generating function associated the sequence (a-,), we first note that the sequence 
(flj) satisfies the recurrence relation 



(3.10) 



if j < 
if J = 0, 
otherwise. 



1 



Qjj \ ~\~ dj pT7 

This recurrence relation can be expressed as the single equation 



(3.11) 



aj aj-i + aj-p,^ + [j = 0] 
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where [j = 0] adds 1 when j = 0. To demonstrate the generating function for A{z), 
we multiply both sides of ()3.1ip by and sum over j. Thus, we find 

oo oo oo oo 

j=0 j=0 j=0 j=0 

oo oo 

CO oo 

(3.12) = z ^ a^z^ + zP"" QjZ^ + 1. 

j=o j=o 

Consequently, the generating function for A{z) is given by 1/(1 — z — z^"). We can 
restate our problem to find a closed form for A{z), and thus evaluate [z^]A(z). We 
use the following result from Graham et al. (1994). 

Lemma 3.7. (Graham et al. (199^)) f^^_^^y+i = Yln=o Ct")""^"- 

We seek [z^]yl(z) = [z^]l/(l — z — z^ ). Consider the series, 

oo / ■ I • 

(3.13) l/(l-az)^+i=^^^ + ^ 

and further consider a finite sum of such series 

(3.14) S{z) = + . . . + 

^ ' ^ ' il-a,zy'+' (l-ajvz)*"+' 

Then [z-']S'(z) is the finite sum of coefficients given by 

(3.15) [z3]s{z) = A (^'1 + ^y^+.-.+p^ (^^+ ^y^. 

Let A{z) = -gp^ where B{z) = 1 — z — z^ . We now show that J5(z) has distinct 
roots. Let us suppose that B{z) has a set of repeated roots. Then, we have it 
that B{z) and B'{z) share a set of common roots. Since B'{z) = —1 — p™(z)P'"~^, 
it follows that a repeated root a satisfies B'{a) = —1 — p™'{a)P"^ = 0, and 
consequently, ~^ = — 1/p™. Whence, = —a/p"\ Furthermore, as B{a) 
vanishes, we have it that 1 — a + a/p™ = 0, or equivalently, 1 — a(l — l/p™) = 
from which we deduce a = p'"/(p™ — 1) to be the only candidate for a repeating 
root. Therefore, a = p"^/{p"^ — 1) should satisfy the equation a^"""^ = — 1/p™, as 

a consequence of being a root of B'{z). That is (■pm_\-)p'"-i = Now {p"^)^'" = 

(mod p) while — (p'" — 1)^ ^^ = —1 (mod p) implies that a — p"^/{p"^ — 1) is not 
a root of B'{z). Therefore, B{z) has distinct roots. 

Now writing B(z) in the form (z — 6i)...(z — &pm), and taking reciprocals ak 
of k = l,...,p™, we establish a correspondence with the polynomial (1 — 
aiz) ... (1 — Q!pmz). Thus, A{z) = 1/((1 — aiz) ... (1 — apmz)) may take the form 
/3i/(l — aiz) + • • • + /3pm /(I — Qfpmz), for some f3i, 1 = 1,... Note 1/(1 — az) 
is a special case of Lemma |3 . 71 with i = 0. 

Theorem 3.8. We claim that 

(3.16) W]A{z) = I^Aa^ 

1=1 
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where (3i = —ai/B'{l/ai). 
Proof: 

(3.17) lim (z - l/a/)A(z) = lim (z - l/a/)S(z), 

2-^1/q, z^1/q, 

where S(z) is the special case of equation 1 3 . 1 5 1 wit h ii = 0. Now, it follows that 
lim {z-\/ai)A{z) = lim (z-l/a;)— 3— 

z^l/ai z^l/ai B{Z} 

z -l/ai 



(3.18) 
and. 



lim , „, , , 

z^i/c, B{z)-B{l/ai) 

1 



B'{llai) 



(3.19) 



lim (z-l/ai)S(z) = lim (z-l/a;)y^7- r 

z^l/ai z^l/ai '^(1-afez) 

/?/(z-1/q;,) 

= ^™ 7 71 — 7 

z^i/a, -ai{z - l/ai) 

A. 
-ai' 



Consequently, we deduce f3i = -^Tf-fj^, for Z = 1, . . . since 
(3.20) lim (z-l/ai) 

z^l/a, (1 - ttfcZ) 

vanishes for k ^ I and the result follows. 

We now test our closed form, equation (|3.16p . against the sequence of integers 
arising for instances of the binomial summation aj ~ X]i=o {''~^^ i ~^''*) for cases 
= 4 and = 8. 

Example 3.9. Case = 4. The reciprocals of the roots of 1 — z — z'* are given as 

ai = -.8191725134, 

aa = .219447421 - .9144736630t, 

ag = .219447421 + .9144736630t, 

ai = 1.380277569. 

Since (3i = "^"^ ^^^^ 



(3i = .1305102698, 

h = .1610008758+ .1534011260t, 

/33 = .1610008758- .1534011260i, 

/34 .5474879784. 
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Following result (|3.16p . a closed form for the binomial coefficients Uj is then given 

by 

ttj = (.f305f02698)(-.8f9f725134)^' 

+ (.1610+ .1534i)(.2f94- .9144ty 
+ (.1610 - .1534i)(.2194 + .9144t)^' 

(3.21) + (.5474879784)(1. 380277569)^. 

Maple Input: for j from to 25 do; aj; end; 

Maple Output: 1,1,1,1,2,3,4,5,7,10,14,19,26,36,50,69,95, 

131,181,250,345,476,657,907,1252,1728. 

Example 3.10. Case p"' = 8. The reciprocals of the roots of 1 — z — are given as 
ai = -.9115923535, 
aa = -.6157823065- .6871957511t, 
as = -.6157823065+ .687195751U, 
ai = .1033089835- .9564836042^ 
as = .1033089835+ .9564836042^, 
as = .8522421840 - .6352622030(,, 
ar = .8522421840+ .6352622030t, 

(3.22) as = 1.232054631. 
The set /3/, ^ = 1,. . . ,8, is 

Pi = .06378010282, 

P2 = .06449005934+ .02789285455/,, 

P3 = .06449005934 - .02789285455/., 

Pi = .06911712233+ .06926484155/., 

P5 = .06911712233- .06926484155/., 

Pe = .1188399306+ .1719523210/., 

Pr = .1188399306- .1719523210/., 

P^ = .4313256714, 
Again a closed form for the binomial coefficients aj is given by 
aj = (.06378010282)(-. 9115923535)^' 

+ (.0644+ .0278t)(-. 61578230 - .68719575i)^' 
+ (.0644- .0278t)(-. 61578230 + .68719575i)^ 
+ (.0691 + .0692t)(. 10330898 - .95648360t)^' 
+ (.0691 - .0692t)(.10330898+ .956483600^' 
+ (.1188 + .1719t)(.85224218 - .63526220t)J' 
+ (.1188 - .17190(.85224218 + .63526220t)J' 

(3.23) + (.4313256714)(1.232054631)^. 

Maple Input: for j from to 25 do; aj] end; 

Maple Output: l, 1,1, 1,1, 1,1, 1,2,3,4,5,6,7,8,9, 11, 14,18, 23, 29, 

36,44,53,64,78. 
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Having found a closed form for A{z) = '^j-^Q ajz^ , we now outline the classi- 
cal means by which the cycle length of (aj) mod p™ may be determined. Fix 
a particular value of to > 1 and having obtained the closed form for the coeffi- 
cients of the corresponding polynomial A{z), evaluate the closed form of aj for 
j = 0, . . . ,p'"-i(p2™ _ 2). To show that p™-i(p2m _ ^j^g ^^^^^ length of 

(flj) mod p™, show that firstly apm-i(^p2m_i^ = 1 and apm-i(j,2m_i)_j vanishes for 
j ~ 1, ... — 1, and lastly, show the subsequence of p'" opening Is does not 
appear until j is p»"-i(p2m _ Table [1] gives the cycle length of the sequence (aj) 
mod p™ for small values of p™. Cycle lengths for all prime powers up to 3125 have 
been confirmed and they all agree with the conjecture. 

We now consider the sequence (a,) mod d where ^ Y.\=o for 
composite d with prime factorization d = p™^ . . . p™'' . We calculate the cycle 
length of {aj) modulo d by calculating its cycle length modulo p™* for each t = 

1, . . . ,r. The cycle length of X]i=o '"^^'^ pT* ^ foi' ^ ~ 1, ■ • ■ is given 

by the j for which J2i=i ("'"'l"^''') niod p™' equals 1 and which has a preced- 
ing sequence {aj^d+i, ■ • • , o-j-i) = (0, . . . , 0) mod p™*. Given that the mapping 
^dXp"'^ ■■■p"''') ■ ^ ^p'"i X ••• X Zp^r is well-defined then the cycle length 
X]i=o ("'^^'i^^'*') iTiod d of the recurrence relation of (aj) mod d is given by the 
LCm| J2i=o mod pT' y ^- Table □ gives some initial values for the 

cycle length of J2i=o ("'"'"'i"^''*) mod d. We state this result as a theorem. 
Theorem 3.11. Let Ei=o C^^I^^^O ^^^^ ^ ^e an integer sequence and consider 
the decomposition |X]i=o ("'"'" t"^^') "^'"^ ^'t"' } X]i=o ("'"'' t"^^') "^""^ 



6e f/ie cj/c/e length 



X]i=o ("' '''i ^'''') mod pY""* I &e </ie cj/c/e 



direct product of disjoint cycles. Let 
lengths o/ |E^£g ("'~^'^^~^^')moc!p"''| , respectively. Then, the cycle length of 

e-^r^^o -"'^ Ie^o e-'i-'^o mod p^ \ }[^^ . 

The results provided in Table[T]are the cycle lengths of {aj) mod d for some small 
values of d. The results are two-fold; firstly, the results give the minimum number 
of CNOT gates required to effectuate a quantum gate according the periodicity of 
the binomial summation construction. Secondly, use of the binomial summation 
construction illustrates for which dimensions the generalised SWAP gate can be 
realised. In particular, we have a SWAP gate in those dimensions d for which the 
cycle length of the function {oj) mod d is equivalent to —1 mod d. For dimensions d 
where our SWAP gate is not possible then the cycle length of {aj) mod d describes a 
permutation, other than that permutation sought, of the d input qudit states. For 
example, in dimension 4 the cycle length of {aj) mod 4 is equivalent to 2 mod 4. 
Therefore, we have it that the quantum gate associated with this particular instance 
of the binomial summation construction evolves four 4-dimensional quantum states 
such that the state |eo)Q of the first quantum system Aq is transposed with the 
state 162)2 of the third quantum system A2 so that the quantum system ^0 is in 
the state [62)0 and the quantum system A2 is in the state |eo)2- Correspondingly, 
the state \ei)-^ of the second quantum system Ai is transposed with the state 163)3 
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d 








2 


3 = 22 - 1 


3 


8 = 32 - 1 


4 


30 = 2(2* - 1) 


5 


24 = 52 - 1 


6 


LCM(63,728) = 6552 


7 


48 = 72 - 1 


8 


252 = 4(2^ - 1) 


9 


240 = 3(3^ - 1) 



Table 1. Cycle length of J2i=o {^~^\~'^^')mod d. 



of the fourth quantum system ^3 so that the quantum system Ai is in the state 
jes)]^ and the quantum system ^3 is in the state |ei)3. Moreover, the gate network 
for an implementation of a generalised SWAP*-'-* gate, for ^ = 1, . . . , d — 1, is an 
instance of I applications of this design - such a network cycles the initial state 
description through I systems. 

Furthermore, the cycle length of the function aj = X]i=o ("'^f *) i^^od 6 over dimen- 
sion 6 is the least common multiple of the cycle lengths associated with the functions 
ElLl )mod 2 and E -£0 ('t o mod 3. The cycle lengths of E -£0 Ct Vod 2 
and X]i=o (''~/ ')™°d 3 are 63 and 728, respectively. By Theorem 13.111 the cycle 
length of (flj) mod 6 is 6552; the least common multiple of 63 and 728. However, 
since 6552 is equivalent to mod 6, we have it that the CNOT architecture asso- 
ciated with the binomial summation function X]i=o (''"f *)™°d 6 acts trivially on 
the input state of six 6-dimcnsional quantum states. 

4. QUTRIT SWAP 

We improve the WilNOT gate construction (Wilmott and Wild (2008c)) by 
considering a construction based on binomial summations that yields a quantum 
gate which is also composed entirely from uses of the CNOT gate. Fig. [2] gives a 
generalised quantum SWAP gate for qutrit states which is based upon our binomial 
summation construction. By Theorem 2.4 (Wilmott and Wild (20G8c)), let us 
suppose that the first quantum system Aq prepared in the state |eo)o, the second 
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Figure 2. Quantum SWAP gate for qutrit states. 



system Ai prepared in the state |ei)-^ and the third system A2 prepared in the 
state 162)2- Implementing our quantum SWAP gate yields the system ^0 in the 
state |ei)Q, the system Ai is in the state |e2)^ and the system A2 is in the state 
|eo)2 and this is achieved with eight CNOT gates as opposed to the ten CNOT gates 
of the WilNOT construction. Note further that our construction is based on a set of 
binomial coefficients which satisiy the linear recurrence relation whose generating 
function is given by 1/(1 — z — z^* ). Generating functions provide elegant means 
for the storing of information about the nature of coefficients in their expansion 
(Blasiak et al. (2007)). For our purpose, the nature of the recurrence relation 
is seen through periodicity of the coefficients modulo 3. We relate this periodic 
nature to instances of the CNOT gate acting on basis states in determination of 
a quantum swap gate. Furthermore, since our quantum gate is based on a linear 
recurrence relation, we illustrate that knowing part of the quantum gate design will 
completely aid in determining the remainder of the quantum gate just as knowing 
the initial conditions of a recurrence relation determines all other coefficients in the 
recurrence. Such type of construction is unusual in that knowing part of the gate 
does not necessarily yield the remainder. 

We illustrate the action of each use of the generalised CNOT gate arising in our 
construction on the input quantum states and subsequent evolved states of the 
system. The first unitary gate Ui acts on the input state of three arbitrary qutrits 
where the first quantum system Ao on the input state is prepared in the state |eo)g, 
the second system Ai on the input state is prepared in the state |ei)j^ and the third 
system A2 on the input state is prepared in the state 162)2- "^^e action of Ui in the 
input state is given by 



;7i(|a) ® \b) ® \c)) 

= C/i ((ao |0) + ai |1) + a2 |2)) ® {bo |0) + b, |1) + 62 |2)) <8) (cq |0) + c, |1) + C2 |2))) 
= Uiiaoboco |000) + aoboci |001) + aQboC2 |002) + aofoieo |010) + ao^iei |011) 
+ao6ie2 |012) + ao^aeo |020) + ao^sCi |021) + aQb2C2 |022) + aiboco |100) 
+ai6oCi |101) + aiboC2 |102) + ai^ieo |110) + aibiCi |111) + ai6iC2 |112) 
+aib2Co |120) + aib2Ci |121) + 016262 |122) + 02^060 |200) + 036061 |201) 
+026062 |202) + 026160 |210) + 026161 |211) + 026162 |212) + 036260 |220) 
+026261 |221) + 026262 |222)) 
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Figure 3. The unitary matrix Ui 



= aoboco |000) + aobod |001) + aoboC2 |002) + aobico |010) + aobid |011) 
+ao6iC2 |012) + aob2CQ |020) + ao^aCi |021) + ao^aCz |022) + aiboCo |110) 
+ai6oCi |111) + aiboC2 |112) + aibiCo |120) + aibid |121) + aibiC2 |122) 
+aib2Co |100) + ai&2Ci |101) + ai^sCa |102) + asboco |220) + az^oCi |221) 
+a2boC2 |222) + asfoiCo |200) + aafeici |201) + ajfeica |202) + 036200 |210) 
(4.1) +026201 |211) + 026202 |212) . 

The second unitary matrix U2 corresponds to the second CNOT gate given in 
Fig. [2] and acts on the state of the system after apphcation of the first unitary 
transformation. The action of the unitary matrix U2 on the state of the system 

C/i(ao6oCo |000) + ao6oCi |001) + aoboC2 |002) + aobico |010) + aobici |011) 
+ao6iC2 |012) + ao62Co |020) + 006201 |021) + 006202 |022) + ai6oOo |110) 
+ai6oCi |111) + ai6o02 |112) + ai6iOo |120) + ai6iOi |121) + ai6i02 |122) 
+016200 |100) + 016201 |101) + 016202 |102) + 026000 |220) + 026001 |221) 
+0260O2 |222) + 0261O0 |200) + 0261O1 |201) + 0261O2 |202) + 0262O0 |210) 
+0262O1 |211) + 0262O2 |212)) 
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Figure 4. The unitary matrix U2 



U2{aoboCo |000) + ao6oCi |001) + ao5oC2 |002) + aohco |010) + aohci |011) 
+ao6iC2 |012) + aob2CQ |020) + oo^zCi |021) + ao62C2 |022) + aiboCo |110) 
+ai&oCi |111) + ai6oC2 |112) + ai6iCo |120) + a^bici |121) + aibiC2 |122) 
+aib2Co |100) + ai&2Ci |101) + 016202 |102) + a2boCo |220) + 026001 |221} 
+026002 |222) + 026100 |200) + 0261O1 |201) + 026102 |202) + 0262C0 |210} 
+026201 |211) + 026202 |212)) 

= 0060O0 |000) + 0060O1 |001) + 0060O2 |002) + 0061C0 |011) + 0061C1 |012) 
+0061O2 |010) + 0062O0 |022) + 0062O1 |020) + 0062C2 |021) + 0160C0 |111) 
+0160O1 |112) + 0160O2 |110) + 0161O0 |122) + 0161O1 |120) + 0161C2 |121) 
+0162O0 |100) + 0162O1 |101) + 0162O2 |102) + 0260O0 |222) + 0260C1 |220) 
+0260O2 |221) + 0261O0 |200) + 0261O1 |201) + 0261O2 |202) + 0262C0 |211) 
(4.2) +O262Oi|212) + O262O2|210). 

The unitary matrix C/3 is given in Fig. [5] and corresponds to the third CNOT gate 
which acts on the state of the system after appUcation of the second unitary. Its 
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Figure 5. The unitary matrix U3 



action on the state 

U2Ui{aoboCo |000) + oofcoci |001) + ao6oC2 |002) + aohco |010) + ao&iCi |011) 
+aobiC2 |012) + ao&2Co 020) + ao&2Ci |021) + ao62C2 |022) + aiboco |110} 
+ai&oCi |111) + ai&oC2 112) + ai&ico |120) + aibici |121) + aibiC2 |122) 
+aib2Co |100) + aib2Ci 101) + 016202 |102) + 026000 |220) + 026001 |221) 
+0260C2 |222) + 026100 200) + 0261C1 |201) + 026102 |202) + 0262C0 |210) 

+0262Ci|211)+ 0262C2 212)) 

is 

Uaiaoboco |000) + oo6oci |001) + 0060C2 1002) + oo6ico |011) + oo6ici |012) 
+0061C2 |010) + 0062C0 022) + 0062C1 |020) + 0062C2 |021) + oi6oco |111) 
+0160C1 |112) + 0160C2 110) + 0161C0 |122) + 0161 ci |120) + 0161 C2 |121) 
+0162C0 |100) + 0162C1 101) + 0162C2 |102) + 026000 |222) + 0260C1 |220) 
+0260C2 |221) + 0261C0 200) + 0261C1 |201) + 0261C2 |202) + 0262C0 |211) 
+0262C1 |212) + 0262C2 210)) 
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= aoboco |000) + ao^oci |101) + ao6oC2 |202) + aohco |111) + ao^ici |212) 
+aobiC2 |010) + aob2Ca |222) + aob2Ci |020) + ao^aca |121) + ai6oCo |211) 
+aiboCi |012) + ai&oC2 |110) + aibiCo |022) + aibici |120) + ai6iC2 |221) 
+ai62Co |100) + aib2Ci |201) + 016202 |002) + aafeoco |122) + 026001 |220} 
+026002 |021) + 0261C0 |200) + 0261C1 |001) + 0261C2 |102) + 0262C0 |011) 
(4.3) +O262Ci|112) + O262C2|210). 

As outlined, this construction is based on a linear recurrence relation and therefore 
there is a cyclic nature to our gate construction. This is evident when we consider 
our gate within the qutrit setting. As such, our construction is based upon repeat 
application of the unitrary matrices Ui,U2, and U3. As evident in Fig. [2] the fourth 
gate in the construction corresponds to the unitary matrix Ui. In particular, we 
have it that 

l74(?73;72C/i(oo6oCo |000) + Oo6oci |001) + 0060C2 |002) + oo6ico quditOlO + aobici |011) 
+ao6iC2 |012) + 0062C0 |020) + 006201 |021) + 0062C2 |022) + 0160C0 [110) 
+ai6oCi |111) + 0160C2 |112) + O161C0 |120} + 0161C1 |121) + 0161C2 |122) 
+016200 |100) + 0162C1 |101) + 0162C2 |102) + 0260C0 |220) + 0260C1 |221) 
+026002 |222) + 0261C0 |200) + 0261C1 1201} + 0261C2 |202) + 0262C0 |210) 
+0262C1 |211) +0262C2 |212))) 

= C/4(oo6oCo 1000} + 0060C1 |101) + 0060C2 |202) + ao6iCo + ao6iCi |212) 
+0061C2 |010) + 0062C0 |222) + 0062C1 |020) + 0062C2 |121) + 0160C0 [211) 
+0160C1 |012) + 0160C2 |110} + O161C0 |022) + 0161C1 |120} + 0161C2 |221) 
+0162C0 |100) + 0162C1 |201) + 0162C2 |002) + 0260C0 1122} + 0260C1 |220) 
+0260C2 |021) + 0261C0 |200) + 0261C1 |001) + 0261C2 |102) + 0262C0 |011) 
+0262C1 |112) +0262C2 |210)) 

= 0060C0 |000) + 0060C1 |111) + 0060C2 |222) + 0061C0 |121) + 0061C1 |202) 
+ao6iC2 |010) + 0062C0 |212) + 0062C1 |020) + 0062C2 |101) + 0160C0 |201) 
+ai6oCi |012) + 0160C2 1120) + 0161C0 |022) + ai6ici |100) + 0161C2 |211) 
+016200 |110) + 0162C1 |221) + 0162C2 |002) + 0260C0 |102) + 0260C1 |210) 
+0260C2 |021) + 0261C0 |220) + 0261C1 |001) + 0261C2 |112} + 0262C0 |011) 
(4.4) +0262C1 1 122) + 0262C2 |200) . 

In a similar fashion, we have it that 

U5{UiU3U2Ui{aQboCQ |000) + 0060C1 1001) + 0060C2 1002) + oo6ico |010) 
+0061C1 |011) + 0061C2 |012) + 0062C0 |020) + 0062C1 |021) + 0062C2 |022) 
+0160C0 |110) + 0160C1 |111) + 0160C2 |112) + 0161C0 |120) + 0161C1 |121) 
+0161C2 |122) + 0162C0 |100) + 0162C1 |101) + 0162C2 |102) + 0260C0 |220) 
+0260C1 |221) + 0260C2 |222) + 0261C0 |200) + 0261C1 |201) + 0261C2 |202) 
+0262C0 |210) + 0262C1 |211) + 0262C2 |212))) 
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= U5{aoboCQ |000) + aobod |111} + aoboC2 |222) + aofeico |121) + aobia |202) 
+aobiC2 |010) + aob2Co |212) + aobsCi |020) + ao62C2 |101) + ai6oCo |201) 
+aiboCi |012) + ai5oC2 |120) + ai^ico |022) + ai^ici |100) + aibiC2 |211) 
+ai62Co |110) + aib2Ci |221) + aib2C2 |002} + a26oCo |102) + a2boCi |210) 
+a26oC2 |021) + a25iCo |220) + a25ici |001) + a2&iC2 |112) + a2fe2Co |011) 
+a2^>2Ci|122) + 026202 |200)) 

= aofeoco |000) + aoboci |112) + ao6oC2 |221) + ao&iCo |120) + aobici |202) 
+ao6iC2 |011) + aob2Co |210) + ao62Ci |022) + 006202 |101) + a^boCa |201) 
+ai6oCi |010) + ai6oC2 |122) + aibico |021} + aibid |100) + ai6iC2 |212) 
+ai62Co |111) + 016201 |220) + 0162C2 |002) + 0260C0 |102) + 0260C1 |211) 
+026002 |020) + 0261 Co |222) + 0261C1 |001) + 0261C2 |110) + 0262C0 |012) 

(4.5) +026201 1121} + O262O2I2OO) 

where the unitary matrix U5 is given by the matrix U2- The next unitary matrix 
[/g, and as illustrated in Fig. [2 given by this construction is the matrix which 
corresponds to the matrix U3. We then have it that 

Ue{U5UiU3U2Ui{aoboCo |000) + ao6oOi |001) + 006002 |002) + ao6iCo |010) 
+0061C1 |011) + 0061C2 |012) + 0062O0 1020} + 0062C1 |021) + 006202 |022) 
+0160O0 |110) + 0160O1 |111) + 0160O2 |112) + 0161C0 |120) + ai6ici |121) 
+0161 C2 |122) + 0162O0 |100) + 0162O1 |101) + 0162O2 |102) + 026000 |220) 
+026001 |221) + 0260O2 |222) + 0261C0 1200} + 0261C1 |201) + 0261C2 |202) 
+0262O0 |210) + 0262O1 |211) + 026202 |212))) 

= ?76(oo6oOo |000) + 0060O1 |112} + 0,06002 |221) + oo6ico |120) + 0061C1 |202) 
+0061C2 |011) + 0062O0 |210) + 0062O1 |022) + 0062O2 |101) + 0160O0 |201) 
+0160O1 |010) + 0160O2 |122) + 0161C0 1021} + 0161C1 |100) + 0161C2 |212) 
+0162O0 |111) + 0162O1 1220} + 0162O2 1002} + 0260O0 |102) + 0260O1 |211) 
+0260O2 1020} + 0261 Co 1222} + 0261C1 1001} + 0261C2 |110) + 0262O0 |012) 
+0262C1 1121} +026202 |200}) 

= 0060O0 |000} + 0060O1 |012} + 0060O2 |021} + 0061C0 1 120) + 0061C1 1 102) 
+0061C2 |111} + 0062O0 1210} + 0062O1 1222} + 0062O2 |201) + 0160O0 |001) 
+0160O1 1010} + 0160O2 1022} + 0161C0 |121} + 0161C1 |100) + 0161C2 |112) 
+0162O0 1211} + 0162O1 1220} + 0162O2 1202} + 0260O0 |002) + 0260O1 |011) 
+0260O2 1020} + 0261C0 1122} + 0261C1 1101} + 0261C2 |110) + 0262O0 |212) 

(4.6) +0262O1 1221} + 0262O2 1200} . 

Similarly, 

Uj{UQU5UiU3U2Ui{aoboCo |000) + 0060O1 1001} + 0060O2 1002} + 0061C0 |010} 
+0061C1 |011) + 0061C2 |012) + 0062O0 |020) + 0062O1 1021} + 0062O2 1022} 
+0160O0 |110) + 0160O1 |111) + 0160O2 |112) + 0161C0 |120} + 0161C1 |121} 
+0161C2 |122) + 0162O0 |100) + 0162O1 |101) + 0162O2 1102} + 0260O0 1220} 
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+a2boci |221) + 026002 |222) + a2&ico |200) + a2bici |201) + a2biC2 |202) 
+026200 |210) + 026201 |211) + 026202 |212))) 

= ?77(oo6oOo |000) + 0060O1 |012) + 006002 |021) + ao6iCo |120) + 0061O1 |102) 
+006102 |111) + 0062O0 |210) + 0062O1 1222} + 006202 |201) + ai6oCo |001) 
+ai6oOi |010) + 016002 |022) + 0161O0 |121) + 0161O1 |100) + 016103 |112) 
+016200 |211) + 0162O1 |220) + 0162O2 1202} + 0260O0 |002) + 026001 |011) 
+0260O2 |020) + 0261O0 |122) + 0261O1 |101) + 0261O2 |110) + 026200 |212) 
+026201 |221) + O262O2|200)) 

= 0060O0 |000) + 0060O1 |012) + 0060O2 1021} + 0061O0 1 100) + 0061O1 1 112) 
+0061O2 |121) + 0062O0 |200) + 0062O1 |212) + 0062O2 |221) + 0160O0 |001) 
+016001 |010) + 0160O2 |022) + 0161O0 |101) + 0161O1 |110) + 0161O2 |122) 
+0162O0 |201) + 0162O1 |210) + 0162O2 |222) + 0260O0 |002) + 0260O1 |011) 
+0260O2 |020) + 0261O0 |102) + 0261O1 |111) + 0261O2 |120) + 0262O0 |202) 

(4.7) +0262O1 |211) +0262O2 |220) 

where the unitary matrix Ut is the matrix given by Ui. Our final gate as outhncd 
in this construction and evident in Fig. [2] is the unitary matrix [/§. The unitary 
matrix t/g corresponds to the unitary matrix U2, and we have 

UsiU7U6U5UiU3U2Ui{aoboCo |000) + 0060O1 |001) + 0060O2 |002) + 0061O0 |010) 
+0061O1 |011) + 0061O2 |012) + 0062O0 |020) + 0062O1 |021) + 0062O2 |022) 
+0160O0 |110) + 0160O1 |111) + 0160O2 |112) + 0161O0 |120) + 0161O1 |121) 
+0161O2 |122) + 0162O0 |100) + 0162O1 |101) + 0162O2 |102) + 0260O0 |220) 
+0260O1 |221) + 0260O2 |222) + 0261O0 |200) + 0261O1 |201) + 0261O2 |202) 
+0262O0 |210) + 0262O1 |211) + 0262O2 |212))) 

= ?78(oo6oOo |000) + 0060O1 1012) + 0060O2 |021) + 0061O0 |100) + 006101 |112) 
+0061O2 |121) + 0062O0 |200) + 006201 |212) + 0062O2 |221) + 0160O0 |001) 
+0160O1 |010) + 0160O2 |022) + 0161O0 |101) + 0161O1 |110) + 0161O2 |122) 
+0162O0 |201) + 0162O1 |210) + 0162O2 |222) + 0260O0 |002) + 0260O1 |011) 
+0260O2 |020) + 0261O0 |102) + 0261O1 |111) + 0261O2 |120) + 0262O0 |202) 
+0262O1 |211) +0262O2 |220)) 

= 0060O0 |000) + 0060O1 1010) + 0060O2 1020) + 0061O0 1 100) + 0061O1 1 110) 
+0061O2 |120) + 0062O0 |200) + 0062O1 |210) + 0062O2 |220) + 0160O0 |001) 
+0160O1 |011) + 0160O2 |021) + 0161O0 |101) + 0161O1 |111) + 0161O2 |121) 
+0162O0 |201) + 0162O1 |211) + 0162O2 |221) + 0260O0 |002) + 0260O1 |012) 
+0260O2 |022) + 0261O0 |102) + 0261O1 |112) + 0261O2 |122) + 0262O0 |202) 

(4.8) +O2620i|212)+O26202|222). 

Now, the state of the system after apphcation of the final unitary matrix Us may 
be written 

0060O0 |000) + 0060O1 |010) + 0060O2 |020) + 0061O0 1 100) + 0061O1 1 110) 
+0061O2 |120) + 0062O0 1200) + 0062O1 |210) + 0062O2 1220) + 0160O0 |001) 
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Figure 6. The qutrit SWAP matrix 



+ai6oCi |011) + ai&oC2 |021) + aibiCo |101) + aibid |111) + aibiC2 |121) 
+aib2Co |201) + ai&2Ci |211) + ai^sCa |221) + aa^oco |002) + az^oCi |012) 
+a2boC2 |022) + aa&iCo |102) + aafcici |112) + aa^ica |122) + 036200 |202) 
+02^201 |212) + 036202 |222) 
= bocoao |000) + bocoai |001) + 600003 |002) + 6oCiao |010) + 6oCiai |011) 
+600102 |012) + 600200 |020) + 6002O1 |021) + 60C2O2 |022) + 6iCoao |100) 
+6100O1 |101) + 6100O2 |102) + 6101O0 |110) + 61C101 |111) + 610102 |112) 
+610200 |121) + 610201 |121) + 610202 |122) + 62C0O0 |200) + 62C001 |201) 
+6200O2 |202) + 6201O0 |210) + 6201O1 |211) + 62C102 |212) + 62C2O0 |220) 

(4.9) +6202O1 |221) + 6202O2 |222) . 

In particular, we note that state 14.9] has the form 

(60 |0) + 61 |1) + 62 |2)) ® (co |0) + 01 |1) + 02 |2)) ® (00 |0) + 01 |1) + 02 |2)) 

(4.10) ^ |6) ® |o) (g) |o) , 

thereby iUustrating our binomial summation construction to be a quantum SWAP 
gate. 
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